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1 . INTRODUCTION 

\; 

"*  in  the  theory  of  the  cohesive  energy  of  ionic  or  nearly  ionic  crys- 
tals, the  electrostatic  energy,  or  Madelung  energy,  is  important. 
Generally,  the  electrostatic  cohesive  force  persists,  despite  the  pres- 
ence of  strong  repulsive  forces  of  ions  of  similar  charge  among  the 
constituent  ions.  Because  of  the  magnitude  of  the  Madelung  energy,  it 
dominates  the  energy  of  formation  of  most  ionic  solids  and  therefore  is 
of  interest  to  chemists  who  are  attempting  to  grow  single  crystals 
larger  than  microscopic  powders. 

The  amount  of  literature  on  the  various  methods  of  calculating  the 
electrostatic  energy  of  an  ionic  solid  is  absolutely  overwhelming,  but  a 
methodical  compilation  of  the  results  of  these  calculations  does  not 
seem  to  exist.  Practically  every  textbook  of  physical  chemistry  or 
solid-state  physics  lists  a number  of  references  to  original  papers  on 
this  subject.  Consequently,  we  list  here  only  references  pertinent  to 
the  problem  at  hand:  the  calculation  of  the  Madelung  energy  for  Y2HfSs. 
The  technique  used  here  is  elementary;  however,  it  seems  not  to  be  in 
any  of  the  textbooks. 

t\ 

2 . CRYSTAL  STRUCTURE 

The  crystal  structure  of  Y2HfSs  is,  according  to  Jeitschko  and 
Donohue,1  orthorhombic,  Pnma»  which  is  space  group  No.  62  in  the 
International  Tables  of  Crystallography. 2 As  reported  by  Jeitschko  and 
Donohue,  the  cell  dimensions  and  the  atomic  position  of  each  independent 
constituent  ion  are  given  in  table  I.  The  symmetry  operations  necessary 
to  develop  the  position  of  the  ions  in  a single  cell  are  taken  from  the 
International  Tables  of  Crystallography. 2 The  number  of  molecular 
units,  Z,  in  a cell  is  4 so  that  there  are  32  independent  x,  y,  and  z 
coordinates  necessary  to  specify  the  position  of  all  the  ions  in  a 
single  cell.  The  position  of  each  of  the  atoms  is  shown  in  table  II  in 
a coordinate  system  so  that  the  yttrium  ion  in  position  1 is  chosen  as 
the  origin  and  the  ions  are  listed  in  the  sequence  chosen  by  Jeitschko 
and  Donohue. 1 By  translation,  any  lattice  point  is  equivalent  to  any 
other  point  obtained  from  that  point  by  adding  an  integer  to  any  or  all 
of  its  x,  y,  and  z coordinates.  Using  this  rule,  one  can  choose  a 
coordinate  system  so  that  -1/2  < x or  y or  z 1/2,  and  this  system  has 
been  chosen  in  table  II  (p.  7) . 


V.  Jeitschko  and  P.  C.  Donohue,  Acta  Crystallogr . , B31  (1975),  1890. 
2 International  Tables  of  Crystallography,  J,  Kynach  Press, 

Birmingham,  England  (1952),  151. 
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TABLE  I.  ATOMIC  POSITIONS  OF  CONSTITUENT  IONS  IN  Y2HfS5* 


Ion 

Pos 1 t ion 

X 

y 

z 

Y 

8(d) 

0.1778 

0.9974 

0.0251 

Hf 

Mm) 

0.0060 

1/4 

0.5742 

Si 

8(d) 

0. 1*081 

0.0367 

0.1630 

s2 

Mm) 

0.1822 

1/4 

0.3331 

s3 

Mm) 

0.5032 

1/4 

0.5522 

Sa 

Mm) 

0.2921 

1/4 

0.8152 

•Reported  in  W.  Jeitschko  and  P.  C.  Donohue, 
Acta  Crystallogr . , B31  (1975),  1890. 

Note:  The  x,  y , and  z coordinates  are  in 
units  of0a,  b,  and  c,  respectively,  with0a  * 
11.4585  A,  b = 7.7215  A,  and  c = 7.2207  A. 

The  symmetries  of  the  various  positions  are 
8(d),  ±(x,  y,  z),  ±(1/2  + x,  1/2  - y,  1/2  - z) , 
±(x,  1/2  - y,  z) , ±(1/2  + x,  y,  1/2  - z) ; 4 (m) , 
±(x,  1/4,  z),  ±(1/2  - x,  3/4,  1/2  + z). 


3 . COMPUTATION 

The  electrostatic  potential,  at  a site,  y,  due  to  ions  of  charge, 
q^,  in  the  crystal  is  given  by 


4>(u) 


32 

L L 


2,m,n  v«l 


R„ 

mn,vy 


(1) 


where 


R 


2 

imn,vy 


(*„  - \ * *)V  * (y„  - yv  * ")V  * (\  - \ * ”)2°2  ' 


with  l,  m,  and  n integers  ranging  from  positive  to  negative  values  large 
enough  so  that  the  sum  converges  to  a unique  limit  (~10  for  four 
significant  digits).  The  charge  on  each  ion,  Z v,  is  in  units  of  the 
electronic  value  so  that  - eqv.  For  the  normal  valence, 


SOLID  Y2HfS5 


Dn  No. 

V 

Ion 

Coord i nate 

Charge 

qv 

X 

V 

Vv 

z 

V 

i 

Y 

0.0000 

0.0000 

0.0000 

3 

2 

Y 

0.5000 

-0.4948 

0.4498 

3 

3 

Y 

-0.3556 

-0.5000 

-0.0502 

3 

4 

Y 

0. 1 444 

0.0052 

0.5000 

3 

5 

Y 

-0.3556 

0.0052 

-0.0502 

3 

6 

Y 

0 . 1 444 

-0.5000 

0.5000 

3 

7 

Y 

0.0000 

-0.4948 

0.0000 

3 

8 

Y 

0.5000 

0.0000 

0.4498 

3 

9 

Hf 

-0.1718 

0.2526 

-0.4509 

4 

10 

Hf 

-0.1838 

-0.2474 

0.4007 

4 

1 1 

Hf 

0.3162 

-0.2474 

0.0491 

4 

12 

Hf 

0.3282 

0.2526 

-0.0993 

4 

13 

Si 

0.2303 

0.0393 

0.1379 

-2 

IL 

Si 

-0.2697 

0.4659 

0.3119 

-2 

15 

Si 

0.4141 

-0.4607 

-0. 1881 

-2 

16 

Si 

-0.0859 

-0.0341 

-0.3621 

-2 

17 

S! 

0.4141 

-0.0341 

-0.1881 

-2 

18 

S1 

-0.0859 

-0.4607 

-0.3621 

-2 

19 

Si 

0.2303 

0.4659 

0.1379 

-2 

20 

Si 

-0.2697 

0.0393 

0.3119 

-2 

21 

S2 

0. 0044 

0.2526 

0.3080 

-2 

22 

S2 

-0.3600 

-0.2474 

-0.3582 

-2 

23 

S2 

0.1400 

-0.2474 

-0.1920 

-2 

24 

S2 

-0.4956 

0.2526 

0.1418 

-2 

25 

s3 

0.3254 

0.2526 

-0.4729 

-2 

26 

S3 

0.3190 

-0.2474 

0.4227 

-2 

27 

s3 

-0.1810 

-0.2474 

0.0271 

-2 

28 

S3 

-0.1746 

0.2526 

-0.0773 

-2 

29 

S4 

0.1143 

0.2526 

-0.2126 

-2 

30 

s4 

-0.4699 

-0.2474 

0.1624 

-2 

31 

S4 

0.0301 

-0.2474 

0.2874 

-2 

32 

S4 

-0.3587 

0.2526 

-0.3376 

-2 

[Vie  x,  y,  and  z coordinates  are  in  units  of  a, 

■>,  and  c so  that  the  actual  distance  from  ion 

1 to  any  ion,  v,  is  ( a2x 2 + b2y 2 + c2z2^. 

V v v/ 
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q^  = 3,  1 < V < 8 (Y3+) , 

= 4,  9 _<  v < 12  (Hf4+)  , 

qv  = -2  13  < v <_  32  (S2")  , 

32 

and  qv  = 0 for  neutrality  of  the  entire  solid.  The  sum  given  by 
equation  (1)  was  calculated  for  each  of  the  different  ions  (Y,  Hf,  Si, 
S2#  S3,  and  S4)  given  in  table  I.  The  sum  was  done  so  that,  for  each 
value  of  l,  m,  and  n,  the  sum  over  all  the  constituent  ions  in  a cell 
was  done  first.  The  results  of  this  calculation  are  shown  in  table  III. 
As  a check  on  the  convergence  of  the  sum  for  a given  site,  V*  was 
computed  and  compared  to  a site  that  is  identical,  that  is,  4>(1) 


TABLE 

III. 

ELECTRIC  POTENTIAL  AT  VARIOUS 
Y2HfS5;  UNITS  ARE  (A)'1 

SITES  IN 

Ion 

Type 

u 

<t>(u) . qs  = -2 

0 (u) . qs  = 0 

♦ l 

Y 

1 

-2.2256 

-1.7301* 

0.21*760 

Y 

5 

-2.2256 

-1.7302 

0.21*770 

Hf 

9 

-2.1*535 

2.9775 

2.71555 

Si 

13 

1 . **336 

-0.72753 

-1.08057 

Si 

17 

1 .1*337 

-0.72777 

-1 .08071* 

s2 

21 

1 .A088 

0.13992 

-0.631*  1*1* 

S3 

25 

1.U92 

1 . 1821* 

0.03160 

S4 

29 

1 . 600A 

0.11696 

-0.71*172 

Note:  The  charge  on  yttrium  ion  qy  is  +3,  while 

the  charge  on  the  hafnium  ion  is  q^f  chosen 
so  that  qHf  + 5qs  = -6,  where  qs  is  the 
charge  on  the  sulfur  ion.  The  last  column 
is  the  constant  in  the  expression  = 

1 • 
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and  $(2)  or  <P  (9)  and  <J>  ( 10 ) , and  the  sum  was  extended  until  the  agreement 
was  satisfactory.  A limit  of  ±10  on  the  Z,  m,  and  n sums  gave  agreement 
to  four  digits.  The  total  electrostatic  energy,  UE,  for  a cell  is  given 
by 


u 


E 


6 

E <3.N.  4>  (p  - ) , 

i=l  1 1 1 


(2) 


where  is  the  number  of  each  type  of  ion  (such  as  q^  = 3,  =8  for 
Y and  q^  = 4,  = 4 for  Hf)  and  y covers  the  distinct  types  of  ions. 

In  the  discussion  of  the  structure,  Jeitschko  and  Donohue1  state 
that  the  compound  YzHfSs  is  semi-ionic.  That  is,  the  yttrium  site  is 
ionic  along  with  the  complex  (HfSs)6  , but  within  this  latter  complex 
the  hafnium  and  sulfur  may  be  covalent.  To  account  partially  for  this 
covalency,  one  can  choose  an  effective  charge  on  the  sulfur  ion,  qs,  so 
that 


qHf  + 5qs  = ~6  ' (3) 

which  retains  the  ionic  character  of  the  HfSs  complex.  Also,  the 
lattice  sum  was  performed  by  using  qs  = -0  (qHf  = -6,  not  realistic), 
and  the  results  are  given  in  table  III.  Since  the  electrostatic  poten- 
tial is  linear  in  the  charge,  each  potential  can  be  written  in  the  form 


<p  (y)  = 4>o  (u)  + qg<h  (y) 


(4) 


for  each  distinct  site.  In  the  simple  compounds  where  the  x-ray  data 
depend  on  only  one  parameter  such  as  "a,"  the  cell  size  in  a simple 
cubic  material,  the  electrostatic  potential  is  generally  expressed  as 


$ = 


(5) 


where  M is  the  Madelung  constant,  which  has  been  tabulated  for  most  of 
the  simpler  cubic  structures.3  In  the  more  complicated  structures,  an 
effective  Madelung  constant  can  be  defined  by 


V.  Jeitschko  and  P.  C.  Donohue,  Acta  Crystallogr. , B31  (1975),  1890. 
3 Mendel  Sachs,  Solid  State  Theory,  McGraw-Hill  Book  Co.,  New  York 
(1963),  55. 
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= ((i(y)v1 */3  , 


(6) 


where  <J>(p)  is  given  by  Jeitschko  and  Donohue,1  and  v is  the  volume  of 
the  unit  cell.  We  do  not  use  this  definition  here. 

If  one  chooses  a solid  of  dimension  A,  B,  and  C lengths  in  the  x,  y, 
and  z dimensions,  respectively,  then  the  number  of  cells  in  its  volume 
is  given  by 


N 


(7) 


and  the  total  energy  per  unit  volume,  w£,  is 


(8) 


where  we  have  used  equation  (2)  and  v = abc. 

In  many  calculations,  where  the  complete  binding  energy  is  required, 
one  must  assume  some  form  of  a repulsive  force.  Many  different  types 
cure  assumed  in  the  literature  with  various  degrees  of  sophistication. 
We  assume  the  simple  Born  potential,  B/rn,  where  B is  a constant  and  n 
is  determined  from  experimental  data.4  If  the  material  is  assumed 
isotropic,  then  both  the  electrostatic  energy  and  the  repulsive  energy 
can  be  written  in  terms  of  a single  variable  such  as  v1'3,  where  v is 
the  volume  of  the  unit  cell.  The  condition  that  the  energy  is  a minimum 
then  requires  that  the  derivative  of  the  energy  vanish  at  the  equilib- 
rium value  of  v.  If  it  does,  then 


(9) 


with  w given  by  equation  (8) . 

E 

It  is  customary  to  express  the  energy  in  kilocalories  per  mole,  WE; 
for  Y2HfS5,  this  conversion  of  the  various  constants  gives 


1W.  Jeitschko  and  P.  C.  Donohue,  Acta  Crystallogr. , B31  (1975),  1690. 

4A.  J.  Dekker,  Solid  State  Physics,  Prentice-Hall , Inc.,  Englewood 

Cliffs,  NJ  (1965),  Ch  5. 
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WE  = 83.05we 


(10) 


if  wE  is  in  reciprocal  angstroms.  Using  the  results  given  in  table  III, 
one  can  write  electrostatic  energy  as 


w (q  ) = -56.495  - 59.420q  - 34.167q2  . 

E s s s 


(11) 


The  result  given  in  equation  (11)  for  several  values  of  qs  is  wE(-2)  = 
-74.325,  wE(-l)  = -32.240,  and  w£  = -40.275,  all  in  reciprocal 
angstroms.  With  the  result  of  equation  (11)  substituted  in 
equations  (10)  and  (9),  one  finally  has  for  W,  the  total  energy. 


W = 83.05(1  - 


Generally,  the  best  fit  value  of  n is  obtained  by  matching  the  theoret- 
ical and  experimental  values,  and  in  simple  compounds  the  resulting 
value  of  n varies  from  6 to  12.  Choosing  a value  of  n = 10  in 
equation  (12) , one  obtains 


V(- 


56.495  - 59.420q  - 34.167q2)  kcal/mole.  (12) 


W = -4222.8  - 4441. 4qg  - 2553. 8q2  kcal/mole  . (13) 

This  gives  the  lattice  energy  frequently  used  in  the  Born-Haber  cycle 
calculations. 5 

Occasionally,  the  lattice  energy  is  calculated  by  assuming  that  the 

solid  is  formed  from  preexisting  ionic  complexes.6  A reasonable 

approach  from  this  assumption  is  2Y3*  + (HfSs)6",  which  forms  the  solid. 

In  this  approach,  one  needs  the  energy  of  formation,  U,  of  a single 

(Hf S5 ) 6 ' . This  energy  as  a function  of  q is 

s 

u[(HfS5)6']  = -46.709qs  - 28. 332q2 (a) _1  . (14) 


SF.  A.  Cotton,  Advanced  Inorganic  Chemistry,  Interscience  Publishers, 
New  York  (1972),  61. 

6A.  B.  Blake  and  F.  A.  Cotton,  Inorganic  Chemistry,  2_  (1963),  906. 
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If  this  result  is  subtracted  from  equation  (11)  and  the  result  is 
multiplied  by  83.05 [1  - (1/n) ] , then  the  energy  of  formation  of  the 
solid  ^2  HfS5  from  the  constituent  ionic  complexes  2Y3+  and  (HfS5)6-  is 


U = -4222.8  - 948. 74q  - 436. 14q2  kcal/mole  . (15) 

s s 
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